Time resolved quantum dynamics of double ionization in strong laser fields 
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Quantum calculations of a 1+1-dimensional model for double ionization in strong laser fields are 
used to trace the time evolution from the ground state through ionization and rescattering to the 
two electron escape. The subspace of symmetric escape, a prime characteristic of nonsequential 
double ionization, remains accessible by a judicious choice of 1-d coordinates for the electrons. 
The time resolved ionization fluxes show the onset of single and double ionization, the sequence of 
events during the pulse, and the influences of pulse duration, and reveal the relative importance of 
sequential and non-sequential double ionization, even when ionization takes place during the same 
field cycle. 
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Double ionization in intense laser fields has been chal- 
lenging because of a yield much higher than derived from 
independent electron calculations, thus demonstrating 
the significance of electron interations (see Q and ref- 
erences therein). High resolution experiments revealed 
that the two outgoing electrons preferably leave the atom 
side by side, with the same parallel momenta 0,0. The 
theoretical understanding and interpretation of this pro- 
cess is still far from being complete. The most accu- 
rate representations of the process, i.e. the exact so- 
lution of the time-dependent Schrordinger equation for 
two electrons in a laser field |3( or S-matrix calculations 
Q are computationally demanding and still do not fully 
represent the experiments. Low-dimensional models fre- 
quently sacrifice the experimentally dominant subspace 
of symmetric escape by restricting the electrons to move 
along a common line (aligned-electron models) EJ, or 
introduce other correlations, as in 7], where the motion 
of the electron center of mass is restricted to be along 
the field polarization axis. The 1+1-dimensional model 
we present here removes these drawbacks, and allows for 
efficient calculations which give time- and momentum- 
resolved insights into the dynamics of the process, from 
the turn on of the field to the final escape of the elec- 
trons. The calculations reveal the time order of events 
and clearly show the significance of the saddle in the 
symmetric subspace for the ionization yield, thus lending 
further support to the ideas derived from studies of the 
classical models . 

The model is motivated by the rescattering scenario, 
While most electrons leave the atom directly and 
contribute to the single ionization channel, some have 
their paths reversed by the field and return to the core. 
The acceleration by the field brings in enough energy so 
that when this energy is shared with another electron 
close to the nucleus, each has enough energy to ionize. 
During the collision with the other electron, a short lived 
compound state is formed which then decays into differ- 
ent possible channels: double ionization, single ionization 
or a repetition of the rescattering cycle. Starting from 



this intermediate state a classical analysis easily yields 
possible pathways to ionization 8]. The classical model 
of non-sequential double ionization (NSDI) suggests that 
the electrons may escape simultaneously if they pass suf- 
ficiently close to a saddle in the symmetric subspace that 
form in the presence of the electric field. As the field 
phase changes, the saddles for this correlated electron 
escape move along lines that keep a constant angle with 
respect to the polarization axis. 

The observation that the saddles move along lines 
through the origin suggests a model where each electron 
is confined to move along this reaction coordinate [Toj . 
The main difference of this model with the aligned elec- 
tron models is that a symmetric motion of the electrons 
is possible and not blocked by electron repulsion. The 
present model is hence able to reproduce tunneling and 
rescattering processes, single and sequential double ion- 
izations and it correctly mimics the correlated electron 
escape. Moreover, because of the restriction to 1+1 de- 
grees of freedom it can be integrated by standard meth- 
ods. 

Taking into account that the lines form an angle of 7r/6 
with the field axis 8] , the restricted classical Hamiltonian 
for the two electrons in the linearly polarized laser field 
is given by (in atomic units) [Toj : 
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where Ti are the electron coordinates along the saddle 
lines. The electric field is F(t) — Ff(t)sm(ut + </>), 
with amplitude F, envelope /(£), frequency uo (here: uo = 
0.06 a.u.) and phase </>. For a fixed time and field F(t) 
the potential energy JIJ has a saddle located in the invari- 
ant symmetric subspace r\ = r2 and pi = P2 at |n| = 

]p\ = Z^/VWW of energy V s (t) = -3 3 / A 2^\F{t)\ 
[ll[. If the electrons pass close to this saddle and suffi- 
ciently close to the symmetric subspace, they can leave 
the atom simultaneously. The saddle defines the bot- 
tleneck for simultaneous escape. Once the electrons are 
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FIG. 1: Geometry of the model. The left frame shows the 
lines along which the electrons can move. The right frame 
shows the division of the configuration space r\-ri into do- 
mains assigned to the neutral atom A, singly charged ions Si 
and doubly charged ions, D*, each with indices 1, . . . , 4. The 
polarization axis points along z. 



outside the barrier they are accelerated by the field and 
escape. Any asymmetry of electron motion around the 
saddle can be amplified by the field. Thus, even if they 
escape simultaneously their final momenta can be quite 
different (see Ref. [l(| for examples of classical trajec- 
tories in a static field). In this letter we analyze the 
quantum dynamics of this model by solving the time- 
dependent Schrodinger equation numerically by a Fourier 
methods. The Coulomb singularities in the potential in 
JIJ are smoothed by the substitution 1/x — > 1/Vx 2 + e 
with e = 0.6, which leads to a ground state energy of the 
unperturbed atom of E g = —2.83. 

For the analysis of the outgoing electrons, we follow 
Ref. and define regions in the configuration space that 
correspond to the neutral atom (A), the singly ionized 
(Si) and the doubly ionized atom (Di) (see Fig.[Q. These 
definitions are suggested by practical considerations and 
correspond effectively to a truncation of the long range 
effects of the Coulomb potential at these distances. The 
regions allow us to distinguish between the sequential and 
the non-sequential (simultaneous) double ionization by 
calculating the probability fluxes between the appropri- 
ate regions: i) The population of the singly ionized states 
(SI) at a time t is obtained from the time integration of 
the fluxes from A to Si minus the fluxes from Si to D^; 
ii) The population of NSDI states is obtained from the 
time integration of the fluxes from A to Di and from A 
to D3. iii) Integration of the fluxes from to Dj gives a 
measure of sequential double ionization (SDI) processes. 
The fluxes from A to D2 and from A to D4, correspond 
to anti-correlated double ionization: they give negligible 
contributions to the double ionization process and will 
not be considered further here. Note that the definition 
of the fluxes allows us to distinguish two contributions to 
the instantaneous double ionization yield: electrons may 
pass directly from region A to Di, say, or they may first 
cross over to Si and then to Di. The essential difference 
is that in the first case electron interactions are signif- 
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FIG. 2: Yields for single ionization (circles), sequential double 
(triangles) and non- sequential double (squares) ionizations as 
a function of the field amplitude. The data are obtained for 
the initial phase of the field = and the pulse envelope 
f{t) — sin 2 (7r£/T) where pulse duration T equals 5 field cycles. 
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FIG. 3: Probability fluxes (in arbitrary units) as a function of 
time. Panel (a) shows the field strength for F — 0.16. Panel 

(b) shows the flux related to the single ionization and panels 

(c) and (d) show the non-sequential and sequential double 
ionization yields, respectively. The field has an initial phase 
<f> — 0, a duration of 8 cycles, and is switched on and off 
linearly over 2 cycles. 



icant, whereas in the second case the electrons remain 
sufficiently far apart that their interactions are negligi- 
ble. After a rescattering event both processes can hap- 
pen. This suggests that the definition of NSDI should 
be restricted to processes where electron interactions are 
important, i.e. the direct flux from A to Di and D3. In 
the following we will use the term only in this meaning. 

As a first result we show the ionization yields for the 
different subspaces in Fig. [21 They are calculated from 
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FIG. 4: Time resolved electron momentum distributions cor- 
responding to the field parameters in Fig. El The times equal 
the extrema of the field strengths, indicated by the dashed 
lines in Fig. |3| starting with the one at t = 3.75 cycles (top 
left) and ending with 7.25 (bottom right). 



the fluxes as described in the preceeding section. For 
intermediate fields the NSDI and SDI signals are about 
equal, but for higher fields SDI rises sharply, forming 
the well known knee: in hindsight, it is clear that only 
the sequential double ionization can show the knee, as it 
derives from the strong increase of independent electron 
ionization at high fields. The NSDI-signal, on the other 
hand, seems to saturate for field amplitudes above about 
0.25 a.u. 

The time ordering of the process, resolved into the dif- 
ferent fluxes, is shown in Fig. [3| for a field strength of 
F = 0.16, below the knee. Up to the fourth extremum 
the field is not strong enough to ionize any electrons. 
Shortly thereafter singly charged ions appear, but no 
doubly charged ones. Immediately after the fifth ex- 
tremum a strong single ionization is observed as well as a 
first double ionization signal. Note that the maximum of 
double ionization occurs after the extremal field strength, 
at about the same time as the maximum in the single 
ionization signal. This shows that in contrast to expec- 
tations based on the simple man model p], |^ for the re- 
turning electron, the double ionization events occur when 
the field is still on. 

The momentum distributions of the outgoing electrons 
are obtained following the method proposed in 6] where 
the wavefunction is propagated with all interactions in 
a region of width 400 a.u. x 400 a.u. Outside this do- 
main the wave function is transformed to the momen- 
tum space where time evolution (with neglected Coulomb 
potentials and in the velocity gauge) becomes simply a 
multiplication by a time dependent phase. Fourier trans- 
forms of the parts of the wavefunction in the regions 
|ri | , |r2 1 > 200 a.u. then give the momentum distribu- 
tions in Fig. 0] 

The panels in Fig.^Jare calculated for the same param- 
eters as for Fig. |31and give the momentum distribution at 
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FIG. 5: Rescattering in a 1-d model at F — 0.16. Panel 
(a) shows the return time t r , panel (b) the excess energy at 
the recollision moment t r and panel (c) the tunnelling prob- 
ability, obtained from a semiclassical estimate oc e~ s where 
S is the action of a tunneling trajectory, in unsealed units. 
The abscissa for all panels is the point in time where the 
electrons tunnel through the barried (their initial energy is 
—0.83). Dashed lines in panels (a) and (b) show the results 
of the model with neglected Coulomb potential [J |?J . 



successive extrema of the field. The sequence starts with 
the extremum at 3.75 cycles, as there is no noticable wave 
function amplitude in the range | , |r*2 1 > 200 a.u. for 
earlier extrema. The wave function that gives rise to the 
NSDI near times of t = 2.25 oscillates with the field and 
extends into this space region only about 1.5 cycles later. 
Thus, with the delay taken into account, the first signals 
in the NSDI sector in Fig.0]correspond to the first signals 
in Fig. 03 The momentum distributions in the first four 
panels in Fig. 21 start out very much concentrated along 
the diagonal p\ = P2- This confirms that the bottleneck 
for double ionization are the saddle configurations in the 
symmetric subspace, described in 

After a few cycles, the different ionization signals in 
Fig. [3| and the momentum distributions in Fig. ^experi- 
ence significant spreading and distortion. They then no 
longer reflect the sequence of extrema and rescattering 
events, and the temporal relation between the individual 
processes gets blurred. Moreover, as time goes on, less 
correlated and purely sequential processes become more 
important. This suggests that the structure of the pro- 
cess can best be resolved with short pulses, say up to 3 
field cycles. This is shorter than the pul ses used so farQ, 
but now within experimental reach |l2| . 

The time ordering of the process and in particular the 
presence of the field when the electrons return to the nu- 
cleus in the rescattering event, can also be understood 
from the classical dynamics, as in 9], if the Coulomb 
field is taken into account. To this end we show in Fig. El 
the results from a classical trajectory calculation. As in 
|9| we assume that an electron that tunnels out is re- 
leased with zero momentum at the other side of the po- 
tential barrier. It is then integrated classically until it 
returns to the atom. Since the process involves motion 
of a single electron along the field axis, we can take a 1-d 
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FIG. 6: Final electron momentum distributions for different 
field strengths. The pulse is 8 cycles long with a linear switch 
on and off over 2 cycles. The momentum distributions are 
averaged with a Gaussian of width 0.07 a.u. to model ex- 
perimental resolution (and to remove finite size fluctuations 
from our numerical grid). The field phase is = and the 
amplitudes are F = 0.08 (a), F = 0.16 (b) and F = 0.4 (c). 
In order to show the effects of a varying phase and broader 
smoothing, panel (d) shows results for F — 0.16, but averaged 
over an undetermined phase cj) and a wider Gaussian window 
of 0.2 a.u. The data are the Fourier transforms of the parts 
of the wavefunction (evolved one more cycle after the pulse is 
gone) in the regions |n|, \r2\ > 100 a.u. 

Hamiltonian, H x = p 2 /2 - 1/yV 2 + e + r(V3/2)F sm(ut). 
The electron that tunnels through the Stark barrier starts 
with an energy —0.83 a.u., equal to the energy difference 
between the ground state of a He atom and a He + ion in 
our model 9]. If the field is weak, the electron starts far 
from the core and can aquire considerable energy while 
the field brings it back. However, such processes are very 
unlikely since the tunneling probability is negligible. The 
relevant energy parameter when the electron returns to 
the nucleus is the difference between energy of the two 
electron system E(t r ) and the potential energy of the sad- 
dle V s (t r ) HDIl, defining AE = E(t r )-V s (t r ). The data 



collected in Fig. [5| (corresponding to F = 0.16) clearly 
show that most electrons return while the field is still on 
in agreement with the sequence of events documented in 
Fig. El For smaller F the range of positive excess energy 
shrinks and moves towards larger values of to where the 
tunneling probability is negligible. This indicates that 
as F increases one cannot expect a sharp threshold be- 
havior for the correlated simultaneous escape since the 
contributions from the correlated events grow smoothly 
with field amplitude. 

Fig. shows results for an 8 cycle pulse and different 
field strengths. The figure can be compared with Fig. 1 of 
Ref. obtained in the aligned electron model where due 
to the overestimated Coulomb repulsion the area around 
Pi = P2 is not populated. Here, this region is accessible 
and provides information about the correlated electron 
escape. With increasing F the double ionization signal 
increases, but above the knee (F = 0.4) the strong con- 
tributions in the second and fourth quadrant show the 
strong influence of SDL The distributions for undeter- 
mined phase and wide averaging in Fig. still show the 
strong concentration of the momentum near the diago- 
nal, but the interference structures in the corresponding 
Fig. 03 are washed out. The investigation of these inter- 
ference patterns is the subject of ongoing work. 

The 1+1-dimensional model for double ionization of 
atoms in strong laser fields discussed here gives access 
to the time-dependence of events and the distribution 
of final momenta, including the subspace of symmetric 
escape. The good reflection of experimental observations 
suggests that the model can also be used to study other 
aspects of the ionization process, such as the effects of 
different pulse shapes or wave packet interferences. 
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